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Аннотация 

Введение. Устройства сбора и накопления энергии из внешней среды представляют собой маломощные источ- 
ники электрической энергии, которые активно используются, в том числе в автономных приборах мониторинга 
поврежденного состояния различных конструкций. Рабочим элементом этих устройств является пьезоэлектриче- 
ский генератор (ПЭГ) — преобразователь механической энергии в электрическую. Конструирование ПЭГ свя- 
зано с предварительным построением их математических и компьютерных моделей, с помощью которых произ- 
водится расчет и оптимизация конструкций. Одним из способов моделирования и расчета ПЭГ является разра- 
ботка приближенных методов расчета на основе прикладных теорий. В литературе известны и ранее разработаны 
прикладные теории расчета изгибных колебаний многослойных пьезоактивных пластин. Однако информации об 
изгибно-сдвиговых колебаниях, как инструменте повышения эффективности инженерных расчетов описанных 
конструкций, в научной литературе недостаточно. Целью настоящей работы являлась разработка прикладного 
метода расчета изгибных и сдвиговых колебаний пьезокерамических пластин, в том числе пористых. 
Материалы и методы. В качестве пьезоактивного материала пластины используется пьезокерамика РУТ-4, в 
том числе пористая. При использовании пористой керамики жесткость конструкции уменьшается в большей сте- 
пени, чем пьезомодули, что позволяет получить более эффективный ПЭГ при механическом воздействии. Мате- 
матическая постановка осуществлена в рамках линейной теории электроупругости при поляризации пластины 
по толщине. Боковые стороны пластины электродированы, правая сторона закреплена, а на левой задан гладкий 
контакт в вертикальной стенке. Установившиеся колебания пластины вызываются давлением на лицевые поверх- 
ности пластины или разностью электрических потенциалов на электродах. Для расчета характеристик ПЭГ в ра- 
боте предлагается прикладная теория, основанная на гипотезах о распределении характеристик напряженно-де- 
формированного состояния и электрического поля. 

Результаты исследования. Рассмотрены поперечные колебания пьезокерамической пластины в низкочастотной 
области (ниже первого изгибно-сдвигового резонанса). В силу того, что математическая постановка рассмотрена 
в рамках линейной теории упругости, задача разделилась на сумму двух. В первой учитывалось механическое 
воздействие: на лицевые поверхности пластины действует распределенная нагрузка и поперечная сила на левом 
конце, а потенциалы на электродах равны нулю. Во второй задаче механические нагрузки отсутствовали, но за- 
давалась разность потенциалов на электродах. На основе гипотез о распределении деформаций, механических 
напряжений и электрического потенциала обе задачи были сведены к системе обыкновенных дифференциальных 
уравнений и граничных условий. Сравнение с результатами расчетов методом конечных элементов в пакете 
АСЕГАМ показали адекватность предложенной прикладной теории в низкочастотной области. 
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Соловьев А.Н. и др. Об одном методе расчета изгибных и сдвиговых колебаний пористого пьезоэлемента в низкочастотной области 


Обсуждение и заключение. Поскольку постановка задачи рассматривалась в линейной теории электроупругости и 
изучалась низкочастотная область, в работе удалось задачу об изгибных и сдвиговых колебаниях пластины из по- 
ристой пьезокерамики разделить на две: изгибную — с механическим воздействием при нулевых потенциалах и 
сдвиговую — при задании разности потенциалов и нулевом механическом воздействии. Использованы соответ- 
ствующие гипотезы об изгибе и сдвиге, построены две системы обыкновенных дифференциальных уравнений и 
граничных условий, которые решаются аналитически без использования «тяжелых» конечно-элементных пакетов. 
Для сравнения результатов и подтверждения адекватности предложенного метода проведено конечно-элементное 
моделирование таких задач в специализированном пакете АСЕГАМ. Это сравнение показало, что ошибка в опреде- 
лении смещений и электрического потенциала при использовании этого подхода, в случае задания механических 
нагрузок и разности потенциалов, не превышает 6 %. Разработанный в статье метод может быть применен при про- 
ектировании пьезоэлектрических генераторов накопления энергии в низкочастотной области. 
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епегоу ћаѓ аге асійуеІу иѕеа. Тһе аџёопотоиѕ аеуісеѕ Рог топіќогіпе ће датасе4 сопаійоп оҒуагіоиѕ ѕігисќцгеѕ іпсіойе 
Фет аз е1. Тһе уогкіпо еіетепі оЁ {Безе еуісеѕ 1$ а ріелоеІесігіс репегаіог (РЕС) — а сопуегіег оЁ тесһапіса] 
епегоу шю еІесігіса1 епегоу. Тһе 4ез1еп оЁ РЕС 15 аѕѕосіаќеа уһ ће ргешитагу сопѕігисіоп оЁ Феи таетайса| 
апа сотршег то4е|5, у (ће Бер оЁһісһ Фе саісшабйоп апа орітіғайоп оф ѕігисіцгеѕ 1$ саггіей ош. Опе оЁіће уауѕ 
їо тое апа саісшаѓе РЕС 1$ їо 4еуеор арргохітаѓе саісшайоп теѓіһћоаѕ баѕей оп аррПе4 ћеогіеѕ. Тһе арріеа іћеогіеѕ 
Гог саісшабпе Беп4те угаНоп$ оЁ шиШауег р1ехоасйуе рІаіеѕ аге Кпо\т апа ргеуіоџиѕ1у аеуеІіореа іп ће Шегаге. 
Ножеуег, іп ће ѕсіепіћс Шегаиге Феге 15 пої епоџећ іпѓогтайоп абоџќ епт? апа ѕһеаг у1ібгайопѕ аѕ а 4001 юг 
пиргоуте ће еЁбсіепсу оРепошеегте саісшайопѕ оЁ ће 4езсте4 ѕігисіџгеѕ. Тће објесйуе оЁ 1$ уогК \аз о їеуеІор 
ап арреа те®о4 Гог саісшайпе Беп4ш$ апа ѕһеаг уібгайопѕ оЁ ріегосегатіс рІаѓеѕ, пса те рогоцѕ опез. 
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ућісһ таде 1 роѕ51Ые ю оМаш а тоге ейесібуе РЕС опаег тесһапіса! асйоп. Тһе таетайса| Ғогптшабйоп Уаз 
саггіеа ош ут Фе Яате\могК оЁ ће Ппеаг һеогу оҒ еІесігое1аѕіісіќу ми рІаѓе роІагіхайоп ш ћісКпеѕѕ. Тһе 14е$ оЁ 
е рІаѓіе уеге еІесігойаіѓеа, ће гіеһ $14е уаѕ Яхе4, апа а ѕтооїћ сопќасі іп Фе уегаса] уа! уаѕ ѕеї оп ће 1ей ѕ1ае. 
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роепНа!$ аі ће еІесігойеѕ. То сасиае ће сһагасіегіѕсѕ оР РЕС, ће ашћогѕ ргороѕеа ап аррНе4 ћеогу Базе оп 
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Плеаг Феогу оЁ @азисйу, ће ргоЫет уаѕ йіуіаеа шю Фе зат оѓ бо. Тһе ћгѕ опе юоК шю ассоипі ће тесһапіса1 
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ейесі: а йіѕігібиќеа Іоаа апа а ігапѕуегѕе Гогсе аќ їће Іей епа асіеа оп ће гоп зигЁасез оЁ ће рІаѓе, апа ће роѓепіа1ѕ 
аќ е еІесігойеѕ уеге лего. Іп Фе ѕесопа {а$К, іһеге уеге по тесһапісаІ Іоааѕ, би ће роепйа1 а1егепсе уаѕ ѕеѓ аї Фе 
ејесігодеѕ. Ваѕеа оп һуроеѕеѕ абои Һе ӣїіѕігірийоп ог аеѓогтайопѕ, тесһапіса1 ѕігеѕѕеѕ апа еІесігіс роѓепіа], Бо 
ргоБІетѕ уеге гефисе4 іо а ѕуѕіет оЁ огаіпагу Чегепна! едиайопѕ апа Боипдагу сопаіійопѕ. Сотрагіѕоп УИ Фе 
гези$ оЁсаісшабопѕ Бу ће бпіќе еіетепі тео ш ће АСЕГАМ раскаре зВо\е4 ће адедиасу оће ргорозе4 аррПеа 
Феогу іп ће Іоу-Ёгедиепсу геріоп. 

Ріѕсиѕѕіоп апа Сопсіиѕіоп. Ѕїпсе ће Ғогтшайоп оР Фе ргоБіет уаѕ сопѕійегей іп ће Ппеаг Феогу оҒеІесігоеІаѕіісіќу, 
апа ће Іоу-Ёгедиџепсу геріоп маѕ ѕіџагеа, ће уогк зиссеедеа іп йіу1йіпе Фе ргоЫет оҒ Беп4те-5Веаг угайоп$ оЁ а 
рогоиѕ ріелосегатіс р!аёе шю уо: Беп@ те — у тесһапіса] асйоп ай его роѓепіа1ѕ, апа ѕһеаг — уһеп ѕеііпо (ће 
роепна1 аі егепсе апа лего тесһапіса1 асйоп. Тће соггезроп4те һуроһеѕеѕ абош бепіте апа зВеаг уеге изе4. Тууо 
ѕуѕіегѕ оЁ ог4тагу @1#егепбйа1 едиайопѕ апа Боџпагу соп41юоп$, мһісһ уеге ѕоіуеа апа1уйсаПу уіћоиќ ће изе оЁ 
“һеауу” ћпіќе ејетепќ раскареѕ, уеге сопѕігисѓеа. То сотраге ће гезиН$ апа сопӯгт ће айедџасу оЁ ће ргороѕей 
теѓћоа, ће бйпіќе еетепё тпойеіпо оЁ ѕисһ їаѕКѕ \аз сагпеа ои іп а ѕресіаіхеа АСЕГАМ расКкаре. Тһе сотрагіѕоп 
ѕһоуеа пай ће еггог іп Ӣеѓегтіпіпе аіѕрІасегтепіѕ апа еІесігіс роќѓепіа! уућеп их 15 арргоасћ, ш ће сазе оЁзейте 
тесһапіса1 Іоайѕ апа роепна! 1 #егепсеѕ, 14 пої ехсееа 6%. Тһе тећоа аеуеІореа іп ће рарег сап бе арреа іп ће 
аеѕіоп оЁ ріегоеІесігіс сепегаќогѕ Гог епегоу ѕіогаре іп ће Іом-Ёгедиепсу геріоп. 
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Введение. Пьезоэлектрические генераторы (ПЭГ) используются для преобразования механической энергии в 
электрическую с последующим ее накоплением. Одна из областей применения ПЭГ — создание маломощных 
автономных возобновляемых источников электрической энергии. Рабочим элементом ПЭГ является 
пьезокерамический элемент определенной формы. Форма и тип деформации этого элемента определяют 
пьезомодуль, который характеризует преобразование механической энергии деформации в электрическую. Так 
пьезомодуль д3зз связан с растяжением-сжатием вдоль оси поляризации, з —с такой же деформацией в поперечном 
направлении к этой оси, 415 — со сдвигом. Использование пористой керамики позволяет создавать более 
эффективные ПЭГ. Это связано с тем, что модули упругости пористой керамики с ростом пористости убывают 
значительно сильнее, чем пьезомодули. Таким образом при одной и той же механической нагрузке амплитуда 
деформации у пористой керамики будет больше, следовательно, выходной электрический потенциал также больше. 

Расчет ПЭГ может быть произведен методом конечных элементов, реализованном в пакетах АМЗУ5, 
АСЕГАМ, СОМБОГ и других. Для пьезоэлементов, один или два размера которых значительно меньше других 
(пластины, стержни), могут быть построены прикладные теории расчета на основе гипотез о распределении ме- 
ханического и электрического полей. Без применения «тяжелых» конечно-элементных пакетов прикладные тео- 
рии позволяют моделировать различные устройств на основе пьезоактивных материалов. В качестве таких мате- 
риалов рассматривались пьезоэлектрические, пьезомагнитные и композиционные пьезомагнитоэлектрические. 
Построение этих теорий основано на принятии гипотез о распределении механических, электрических и магнит- 
ных полей. Эти гипотезы связаны с характером колебаний упругих и пьезоактивных элементов ПЭГ. Наиболее 
распространенными конструкциями являются активные и полупассивные биморфы на основе многослойных пла- 
стин, поляризованные по толщине с электродами на лицевых поверхностях, совершающие поперечные изгибные 
колебания. Исследованию устройств со сдвиговой деформацией пьезоэлементов посвящен ряд работ. Электри- 
ческая модель с пьезоэлектрическими определяющими уравнениями режима (15 и модель с одной степенью сво- 
боды объединены для описания характеристик сбора энергии пьезоэлектрического кантилевера в сдвиговом ре- 
жиме в работе [1]. Предлагаемая модель используется для моделирования частотной зависимости выходного пи- 
кового напряжения и мощности. Результаты показывают хорошее совпадение с экспериментом и конечно-эле- 
ментным расчетом в АМЗУ5. В работе [2] разработан пьезоэлектрический преобразователь энергии сдвигового 
режима для использования энергии потока воды под давлением. Он преобразует энергию потока в электрическую 
энергию путем пьезоэлектрического преобразования с колебанием пьезоэлектрической пленки. Разработана мо- 
дель конечных элементов для оценки генерируемого напряжения пьезоэлектрической пленки, которая хорошо 
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согласуется с проведенным натурным экспериментом. Одномерная полностью связанная модель колебания 
балки на основе гипотез типа Тимошенко, которая предоставляет единую общую основу для анализа энергии в 
режиме сдвига и изгиба, представлена в работе [3]. В работе [4] изучается влияние неоднородности свойств пла- 
стины при сдвиговых и крутильных колебаниях ее центральной части. В экспериментальной работе [5] был пред- 
ставлен многослойно-цилиндрический пьезоэлектрический сдвиговый актуатор (МСРЗА), работающий в режиме 
сдвига 415, для прецизионного срабатывания при большой механической нагрузке. Актуатор был изготовлен из 
пьезоэлектрических керамических колец РЬ(7лт,Т1)Оз (РИТ-51), которые были концентрически собраны вместе в 
электрически параллельном соединении с попеременно положительной и отрицательной поляризацией в осевом 
направлении. В работе [6] создан метаматериал из идентичных элементарных ячеек, спроектирован и изготовлен 
искусственный прототип устройства с характерными узорчатыми электродами и расположенными в ряд пьезо- 
керамическими субъединицами, который, как доказано, идеально генерирует синтетическую деформацию сдвига 
грани. При том же напряжении возбуждения наблюдается усиление смещения сдвигового типа более чем на по- 
рядок, по сравнению с предыдущими объемными элементами в режиме 415. В статической постановке в работе [7] 
теоретически установлено поле электромеханической связи в сдвигово-изгибающем режиме для кольцеобразной 
пьезоэлектрической пластины. В соответствии с классической теорией упругих пластин малого изгиба и пьезо- 
электрическими определяющими уравнениями было достигнуто аналитическое решение изгибной деформации 
пьезоактюатора под действием электрического поля и концентрированной или равномерно распределенной ме- 
ханической нагрузки. Механизм создания изгибной деформации объясняется осесимметричной сдвиговой де- 
формацией, которая дополнительно вызывает изгибную деформацию одной пьезоэлектрической пластины в 
форме кольца. Этот механизм существенно отличается от механизма пьезоэлектрических биморфных или уни- 
морфных приводов, о которых сообщалось ранее. Проведена оптимизация конструкции кольцеобразного пьезо- 
актуатора. В работе [8] на основе одномерной модели строится функция отклика датчика на основе сдвиговых 
резонаторов (срезы кварца) объемной акустической волны, которые перспективны для поточных измерений вяз- 
кости жидкости, например, в промышленных процессах. В работе [9] с помощью метода конечных элементов 
исследовался пьезоэлектрический преобразователь управления высотой полета с использованием деформации 
модели сдвига. В [10] теория функционально-градиентной пластины с четырьмя неизвестными сдвиговой дефор- 
мации применяется для выражения компоненты смещения. Распределение электрического потенциала представ- 
ляет собой линейную функцию по толщине. Пластина находится под механической нагрузкой и электрическим 
напряжением. Основные уравнения и граничные условия выводятся с использованием принципа виртуальной 
работы. Проведен анализ напряжений и деформаций от параметров конструкции. Электромеханический анализ 
потери устойчивости пьезоэлектрической нанопластины при сдвиге с использованием модифицированной тео- 
рии парных напряжений с различными граничными условиями изучался в работе [11]. Чтобы учесть электриче- 
ские эффекты, к пьезоэлектрической нанопластине прикладывали внешнее электрическое напряжение. Была ис- 
пользована упрощенная теория сдвиговой деформации первого порядка. Основные дифференциальные уравне- 
ния были получены с использованием принципа Гамильтона и нелинейных деформаций Фон-Кармана. В итоге 
результаты показали, что влияние внешнего электрического напряжения на критическую сдвиговую нагрузку, 
возникающую на пьезоэлектрической нанопластине, незначительно. В работе [12] с помощью комбинации двух 
классических подходов моделирования нелинейного поведения пьезоэлектрических материалов исследуется пье- 
зоэлектрический привод сдвигового типа для атомно-силового микроскопа. В частности, новизна предлагаемого 
метода состоит в том, что он сочетает в себе два источника нелинейности полезависимой модели Мюллера и 
Чжана с частотно-зависимой моделью Дамьяновича. Численные результаты, полученные с помощью метода ко- 
нечных элементов (МКЭ), сравниваются с экспериментом. 

Менее исследованы в научной литературе колебания, в которых кроме изгиба присутствует сдвиг, т.е. «рабо- 
чим» является пьезомодуль 415, значение которого убывает с увеличением пористости, но в меньшей степени, 
чем упругие модули. Последнее обстоятельство позволяет построить эффективное устройство преобразования 
энергии. Поэтому разработка прикладной теории расчета ПЭГ с использованием пористой пьезокерамики на ос- 
нове упрощенных моделей без использования «тяжелых» конечно-элементных пакетов представляется весьма 
актуальной задачей. Целью настоящей работы явилось построение прикладного метода расчета для установив- 
шихся поперечных колебаний в низкочастотной области пористой пьезокерамической пластины, характеризую- 
щиеся как сдвигом, так и изгибом. 

Материалы и методы. Исследуемый ПЭГ представляет собой пьезокерамическую пластину (длина /, 
толщина /), поляризованную по толщине, консольно-закрепленную по правой боковой стороне, левая боковая 
сторона прикреплена к инерционной массе, которая совершает вертикальные колебания и закреплена в 
горизонтальном направлении. Электроды расположены на боковых сторонах пластины, поэтому при разности 
потенциалов на них и отсутствии механической нагрузки, основной деформацией в низкочастотной области 
является сдвиг. Рассматриваются колебания, частота которых меньше частоты первого резонанса. 
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Математическая постановка задачи 
Математическая постановка задачи описывается системой дифференциальных уравнений [13] и соответству- 
ющими граничными условиями. 
Олї + 04р Ус = Ё; У:р=0 
с=с? (#+В42)-е/ -Е Е = -УФф (1) 
р+с̧ад =е,-(&+642)+э7-Е &=(Уи+Уи”)/2 
При рассмотрении пористой керамики связности 3-0 в уравнениях (1) используются эффективные физиче- 


ские константы, определенные с помощью пакета АСЕГАМ-СОМРОЗ [14]. Эти эффективные свойства на основе 
представительных объемов (рис. 1) получены в работе [15] и представлены в таблице 1. 


Рис. 1. Представительные объемы в пакете АСЕГАМ-СОМРОЗ [14] 


Таблица 1 
Эффективные свойства пористой керамики 

% пористости 0 10 20 30 40 50 60 70 80 
р, Ке/т 7500 | 6750 | 6000 | 5250 | 4500 | 3750 | 3000 | 2250 | 1500 
сЕ 1010, МИ? 13,90 | 11,56 | 9,25 6,85 5,05 3,34 2,07 1,26 0,68 
сеї 1010, МИ? 7,78 6,15 4,66 3,14 | 2,10 116 0,62 0,28 0,13 
СЕР, 1010, МИ? 7,43 5,82 4,25 2,82 1,87 1,06 0,52 0,24 0,10 
сеї 1010, №2 11,50 | 9,53 7,23 5,42 3,91 270 1,63 0,91 0,47 
сЕеї 1010, №2 2,56 2,23 1,83 1,44 1,10 0,74 0,44 0,23 0,10 
её", С/т2 15,10 | 13,38 | 11,37 | 9,59 7,68 5,93 3,93 2,30 1,25 
её", Си? -520 | =428. | За | -2 1.2132 | 025 | 504 |0 | 2010 
её", С/т2 12,70 | 10,96 | 8,96 6,91 5,00 3,30 1,95 1,00 | 0,44 
ки Уер 730 663 582 509 439 349 263 191 122 

ке / ео 635 567 492 413 345 270 199 130 75 


По данным таблицы 1 построены зависимости эффективных упругих модулей и пьезомодулей от процента 
пористости [15], которые представлены на рис. 2. 


100 500 
<) > 300 ; 
Е 60 — ЕЁ ЕА ГЕ 200 — ай 
Е — ЕЕ Е — 49 
а 40 з Е 100 а 
а Бо 

%• 4 
си & 100 р 
0 Е 
0 20 40 60 80 в 30 40 50 60 70 80 
РогозНу, % РогозНу, % 
а) 6) 


Рис. 2. Зависимости значений от процента пористости [15]: а — упругих моделей; б — пьезомодулей 
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В соответствии с результатами, представленными на рис. 2, модуль упругости Е! убывает значительно быст- 
рей, чем пьезомодуль 415. 

Построение прикладной теории 
Построенный метод расчета состоит из решения двух задач: в первой рассматривается изгиб под действием ме- 
ханической нагрузки при нулевой разности потенциалов, во второй — моделируется сдвиг, вызванный разностью 
потенциалов при отсутствии механической нагрузки. В обоих случаях учитывается отсутствие зарядов на лице- 
вых поверхностях пластины. При действии механической нагрузки и разности электрических потенциалов ре- 
зультаты двух задач складываются в силу линейности задачи. 

Действие механической нагрузки, потенциалы на электродах равны нулю. Принимаются гипотезы типа 
Кирхгофа-Лява относительно равенства нулю нормальных напряжений, а для перемещений: 


а 
из =02> (х), ш =-| — 022 (х) |2. (2) 
ах 
Распределение электрического потенциала по толщине предполагается квадратичным и симметричным: 
2 2 22 
ф = Ф| (22-1) + Ф (а) 22) + Ф, (х) т Б (3) 
Для учета граничных условий на концах пластины (х = 0, /) получено выражение для поперечной силы: 


8сззозз +8ез 8 


Сізезз 5 7 
2 
1 (сз +ећ) 8сззезз + 8е33 8 13 
31 


О +е 
12 Сзз (сззазз +еђ)л? р? 


1 2 
В 1 8сзз0зз + 8е33 4 р? Р" (20:00) 


Я (сззазз +ез) 12 А? ах 7 
езз (свеззй? — сзземй?) 
С13 д 5 С13 А А 
1 (соза +е%)л ез (сізеззћ — сззезій ) 
+ С + 
12 Сзз (сззғзз +ез) 1? 
2 2 
1 ё5 (свезй — Сззезій )^ а? 
2 7022 (х) 
24 Сзз зз + ез ах 


С учетом равенства нулю нормальной компоненты вектора электрической индукции на лицевых поверхно- 
стях (2 ==), уравнения для неизвестного прогиба 02›(х) и распределения электрического потенциала Ф>(х) 
имеют вид: 


8сззозз +8ез 8 


С1зезз 


(с +е? ) р? 2 
1 33533 33 аа 8сзз 033 + 8е3з 8 13 
31 
2 Сзз (сззазз +ебђ)л? п? 
1 1 8сззезз +8е3 4 |з і 
— — 615 т Эа 92, А —е15й а (х) + 
2 2 (сззёзз +е%)л й ах (5) 
езз (свеззй? сузем?) 
СІЗ > С13 5 з 
1 (сэз +еб)л ез1 (свеззй — сззезіћ ) А 
+ С11 + 
2 Сзз (сззазз +23)? 


1 е15 сзеззй” — сззезаћ? й Г а“ 
| ) 2022(х) – У °рћ022(х)– р(х) = 0, 


24 Сзз зз + езз ах 
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2 8сзз03з + 8е?3 8 
езз ( и Е р? Я 
33533 т 633 8сззезз +8езз 8 
833 5 А 5 Ф› (х) + 
Сзз (сззазз + езз )һ р 


2 2 
+| -2и А ' Е - р? Ф, (х) + 
2 2 (соза ЕЕ еб )л А ах 


6 
(свеззй? = сзземй? ) ( ) 
езз Сіз 
033 (свеззй? —сзземй?) (сззазз +еб)л? а? 
ы 1 7022 (х) - 
(соза че)? Сз3 ах 
а? 
+ Еп (свеззй? — сзземй* Е 8) 
х 
О 2—40), 


5а 2 
24 С338 33 + езз 
Задана разность потенциалов, механическая нагрузка равна нулю. Предполагается независимость попе- 
речного смещения от толщины и равенство нулю продольного перемещения и квадратичное распределение элек- 
трического потенциала по толщине (3): 
Из = (772 (х), ш = 0. (7) 


Выражение для поперечной силы: 
а а 
== —0/ ћ- —Ф р. 8 
О “2 2 (9) (4 2 0) (8) 


С учетом равенства нулю нормальной компоненты вектора электрической индукции на лицевых поверхно- 
стях (2 ==), уравнения для неизвестного прогиба 022(х) и распределения электрического потенциала Ф>(х) 


имеют вид: 
а? а" 2 
=СА4 ——02> (х) й-—е15 —Ф(х) ћ-р(х)- У р^0О2 (х) = 0, (9) 
ах ах 
2 2 
02 (а) ви (а) 0. (10) 


Результаты исследования. Проведено сравнение результатов расчетов по предложенным прикладным 
теориям с расчетами колебаний пьезоэлемента (/ = 0,1 м, л = 0,01 м) на частоте равной 100 с! методом конечных 


элементов в АСЕГАМ [16]. 
В первой задаче, определяемой уравнениями (5) и (6), при задании равномерно распределенной нагрузки 


р(х) = 1000 Па:м и с граничными условиями: 
02 (0) = 0,0 -0= 0, Ф (0) = Ф, (/) = 0, 02 (1) = 0, 502 (7) =0 (11) 


получены следующие результаты, представленные на рис. 3; 4. 
Расчеты показали, что погрешность в определении вертикального смещения — 5,8 %, а для горизонтального 


смещения составляет 1,2 %. 


4,46 1Е-7 Г 
И, м 
27888-7 3х107 
2х107 
1,394Е-7 7 
1х107 
0,0008+0 
0 


0,02 0,06 х, м 
а) 6) 


Рис. 3. Вертикальное смещение в АСЕГАМ в первой задаче: 
а — распределение; 6б — график на верхней границе 
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Рис. 4. Горизонтальное смещение в АСЕГАМ в первой задаче: 
а — распределение; 6б — график на верхней границе 


Во второй задаче, определяемой уравнениями (8) и (9), при задании нулевой нагрузки р(х) = 0, разности по- 
тенциалов Ио = 100 В и граничными условиями: 


О -0= 0, Ф (0) =0,Ф (1) = №, 02 (1) =0 (12) 


получены следующие результаты, представленные на рис. 5, 6. 
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Рис. 5. Вертикальное смещение в АСЕГАМ во второй задаче: 
а — распределение; 6б — график на верхней границе 
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Рис. 6. Электрический потенциал в АСЕГАМ во второй задаче: 
а — распределение; б — зависимость от продольной координаты в середине толщины 


Расчеты показали, что погрешность в определении вертикального смещения составляет 0,8 %, а для электри- 
ческого потенциала — менее 1 %. Следует отметить, что значения горизонтального смещения, рассчитанные в 
АСЕГАМ, оказались на три порядка меньше максимального вертикального смещения, что говорит об адекватно- 
сти гипотезы (7). 

При задании механических нагрузок и разности потенциалов погрешность предложенного метода оказалась 
порядка 6 % для компонент перемещения и электрического потенциала. 
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Обсуждение и заключение. Как отмечалось в цитируемой литературе, одновременное использование изгиба 
и сдвига пьезоэлемента может значительно повысить его эффективность. Кроме этого, использование пористой 
керамики в силу различных зависимостей упругих модулей и пьезомодулей от процента пористости также 
улучшает выходные характеристики ПЭГ. 

В настоящей работе в силу линейной постановки в теории электроупругости удалось построить прикладную тео- 
рию расчета изгибно-сдвиговых колебаний пьезоэлемента в низкочастотной области, которая состоит из решения двух 
задач: в первой действуют механические нагрузки при нулевых потенциалах, а во второй — наоборот — нулевые 
механические нагрузки и задана разность потенциалов. На основе различных гипотез о распределении механического 
и электрического полей получены две краевые задачи для систем обыкновенных дифференциальных уравнений, ко- 
торые решались аналитически. Проведено сравнение результатов расчета смещений и электрического потенциала по 
предложенному методу и с помощью МКЭ, реализованного в пакете АСЕГАМ. Эти расчеты подтвердили примени- 
мость предложенного метода, погрешность для которого в вычислении вышеуказанных характеристик составила 6 %. 
Такой точности достаточно для инженерных расчетов, поэтому предложенный метод может быть применен при про- 
ектировании пьезоэлектрических устройств, в том числе при сборе и накоплении энергии. Дальнейшее развитие этой 
прикладной теории будет направлено на охват более широкого частотного диапазона, включая первый изгибно-сдви- 
говый резонанс. 
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АП аићорѕ һауе геа4 апа арргоуей йе ўпа! тапиѕсгірі. 
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